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1 Introduction. 



More than 20 years have passed since the advent of QCD, but the crucial physical ques- 
tions about the mechanism of confinement and the structure of the QCD vacuum state 
still remain basically unanswered. The QCD sum rule method Q proved to be very useful 
in understanding the hadron spectrum and dynamics, but only rather crude character- 
istics of the vacuum, such as the gluon condensate < G a ^ y G a ^ u > and quark condensate 
< qq >o, could be determined in this way. 

There are two other directions of research which seem now most promising: lattice 
numerical simulations || and model simulations p3] which can be used to confront as- 
sumptions on the form of the QCD vacuum functional with experiment. Naturally, both 
lattice and model simulations are performed in a finite volume. 

Recently, it has been observed that, besides experimental data, there are some exact 
theoretical results specific for finite volume systems which follow from first principles and 
the assumption that chiral symmetry breaking occurs [f|]. In particular, the light quark 
mass dependence of the finite volume QCD partition function has been determined and, 
on basis of that, sum rules for the eigenvalues of Euclidean Dirac operator have been 
derived. The simplest such sum rule is 



where V is the 4-dimensional volume, £ = | < ipip >o |, and v = (1/327T 2 ) / d A xG a ^ LU G a flv is 
the topological charge of the gauge field configuration. The averaging goes over all gauge 
fields with given v with the weight 



where T> is the massless Dirac operator and det' is the product of all its nonzero eigenvalues 
(the index theorem enforces the appearance of v zero fermion modes for each flavor in 
the gluon background with topological charge v. They are shifted from zero when a 



eigenvalues X n which are much less than the characteristic hadron scale oc Aqcd (The 
region A n ^> Aqcd brings about the ultraviolet divergent contribution oc A^ ltr V which 
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(1.1) 




small fermion mass m is switched on). The sum in LHS of ( |1 . 1| ) runs only over positive 



bears no interesting dynamic information). The result ( |1 . 1|) is valid when the length of 
the box L is much larger than Aq 1 cd . 

The sum rule (|1 . 1| ) may be checked by lattice and model simulations, or more correctly, 
be used as a test for the correctness of these simulations. Indeed, if we approximate the 
sum over all gauge configurations by a liquid of instantons, the sum rule (1.1) is satisfied 
1- 

The sum rule (\L.1\) has been derived under the assumption of the standard pattern of 
chiral symmetry breaking 

SU L (N f ) ® SU R (N f ) -> SUv(N f ). (1.2) 

However, the breaking according to this scheme only occurs for fermions belonging to the 
complex fundamental representation of the color group. This is the case when iV c > 3. For 
N c = 2, the fundamental representation is pseudoreal: quarks and antiquarks transform 
in the same way under the action of the gauge group and the pattern of chiral symmetry 
breaking is different leading to different sum rules . A third pattern of chiral symmetry 
breaking is for fermions in the adjoint (real) representation leading to yet another class 
of sum rules |fj. . 

As we shall discuss in more detail in the next section, the true chiral symmetry group 
of the lagrangian of SU(2)— color with fundamental fermions is SU(2Nf) rather than 
just SUi{Nf) (g> SUfi(Nf) (it involves also the transformations that mix quarks with 
antiquarks). The pattern of the spontaneous symmetry breaking due to formation of a 
quark condensate is 

SU(2N f ) - Sp(2N f ). (1.3) 



One of the main results of this paper is the derivation of the analog of sum rule (|1.1|) 

for the case N c = 2: 



= W . (1.4) 

\V^/„ 4(H+2JV / -1) 1 > 

It is of significant practical importance as numerical simulations of QCD with N c = 2 are 
much easier and the sum rule ( p..4|) may be checked sooner. 



In the case of Majorana fermions belonging to adjoint representation of the color 
group has been also considered. The pattern of symmetry breaking is 

SU(Nf) -> SO{N f ). (1.5) 

The analog of the sum rule ( |1 . 1| ) was derived in || for Nf = 1 and Nf = 2 which are 
technically simpler. In this paper we fill up this gap and derive the sum rule for arbitrary 
Nf. 

/ v l\ ^v? nm 

\^KL 4(|P| + (JV / + l)/2)' 1 • ] 

where v = vN c . The sum runs over the positive doubly degenerate pairs of eigenvalues 
of the adjoint Dirac operator. 

Actually, the results ( |1.4|) and ( |1.6|) are not new. They have been obtained earlier 
using universality arguments and random matrix theory ||. In that approach it was 
argued that the spectrum of the Dirac operator near zero virtuality is universal, i.e. it is 
completely determined by the symmetries of the system and can therefore be obtained by 
a random matrix theory with only these symmetries as input. The aim of this paper is to 
illustrate once more that both approaches (stochastic matrices and chiral lagrangian) are 
physically equivalent. The results ( |1 . 1| ) , ( |L4| ) and ( |1.6| ) follow from only the symmetry 
properties of the theory which are properly accounted for in both ways of reasoning. 

The structure of the paper is the following. In the next section, we discuss symmetry 
breaking patterns for QCD with N c = 2 and for the theory with adjoint fermions. In 
section 3, we discuss the technically simpler case N c = Nf = 2, where a nice closed 
expression for the finite volume partition function can be derived. In section 4, we derive 
sum rules ( |1.4|) and (|1.6|) for any Nf and any topological charge sector. In section 5, 
we derive expressions for partition function in the sector with a given topological charge 
for the case of equal fermion masses both for N c = 2 with fundamental fermions and 
for adjoint theories. The results, which generalize the corresponding expression for QCD 
with N c > 3 derived in [[|, can be expressed in terms of a determinant of antisymmetric 
matrices made of Bessel functions (fundamental fermions with N c = 2) or certain integrals 
related to series of Bessel functions (adjoint fermions). 
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2 Chiral symmetry and its breaking. 



In this section we discuss the chiral symmetry breaking pattern for SU(2) color with 
fundamental fermions and for adjoint fermions with SU(N C ), N c > 2. The discussion up 
to (2.9) will be in Minkowsky space. All other sections of this paper and the remainder 
of this section starting from (2.10) will be in Euclidean space time. 

The fermion part of the QCD lagrangian involving Nf massless quark Dirac fields 
belonging to the fundamental representation of the color group is habitually written as 

Cferm = i £ [q f L Vq f L + q f R Vq f R ] , (2.1) 

/=1 

where qL,R = |(1 ± 7 5 )g and D is the Dirac operator. In this basis the U(Nf) <8> U(Nf) 
symmetry of the lagrangian is immediately obvious. Of course, a Ua(1) subgroup is broken 
explicitly by the anomaly. 

However, for N c = 2, the symmetry is higher. To see that, one should note that the 
combination 

~4l = *iC$ (2.2) 

is a left-handed spinor that transforms according to the same representation of the SU(2) 
color group as q( L (i is the color index and C is the charge conjugation matrix). It is 
convenient to write ( [2.1| ) in terms of Nf + Nf two-component Weyl spinors 

2N f 

^ferm = i^W f T>W S . (2.3) 
/=1 

Obviously, this lagrangian enjoys the U(2Nf) symmetry. The U(l) part of it is anomalous 
and only the SU(2Nf) symmetry is left in the full quantum theory. A subgroup is broken 
spontaneously, however, by the formation of the fermion condensate. The condensates 
with maximal flavor symmetry [0] are given by 

(e ij e^w{ a wQ o = |///'. (2.4) 

The matrix I is the 2Nf x 2Nf antisymmetric matrix 

'=(-1 D- m 



where 1 is the Nf x Nf unit matrix. The fact that 1^ is antisymmetric follows just 
from the Grassmann nature of the fields w( a . The specific form ( |2.5|) (It is fixed up to 
an arbitrary transformation / — > UIU T with U G SU(2Nf) , U T is the transpose of 
U. Multiplying of the condensate matrix by an overall U(l) phase factor is also possible 
but it corresponds to going into a sector with a different vacuum angle 9.) follows from 
the requirement that the formation of the condensate does not lead to the spontaneous 
breaking of the vector symmetry SUv{Nf). This is a consequence of the Vafa-Witten 
theorem ||: the spontaneous breaking of a vector flavor symmetry would result in the 
appearance of scalar Goldstone particles. However, the vector-like nature of QCD enforces 
the lightest particle in the spectrum to be a pseudoscalar f\ |5], [H| [□]] . The same conclusion 



follows from an analysis of the effective potential for N c — > oo [O. 

Naturally, the mere fact that the condensate is formed and that the symmetry is broken 
cannot be proven rigorously and is an assumption. For standard QCD with N c = 3 , the 
breaking of chiral symmetry is an experimental fact. It is natural to think that a similar 
(not yet quite disclosed) mechanism also leads to the formation of a chiral condensate 
in the other theories discussed above. The transformations of SU(2Nf) which leave the 
form I?? WfWy invariant constitute the symplectic group Sp(2Nj) (of which SUy(Nf) is 
a subgroup), and we arrive at the chiral breaking pattern ( |1.3|) . 

This result is known in the literature. People were interested with nonstandard pat- 
terns of chiral breaking mainly in association with technicolor models |7], |14| . In a recent 



work ]nj, this question has also been addressed, but the authors assumed a flavor asym- 
metric condensate formation with a different symmetry breaking pattern. 
The breaking (|1.3|) leads to the appearance of 



(2N f y - 1 - (2Nj + N f ) = 2Nj-N f -l 

Goldstone bosons which are parameterized by the coset SU(2Nf)/Sp(2Nf). For Nf = 2, 
we have 5 instead of the usual 3 Goldstone bosons. 

Up to now, we have considered only massless fermions. The mass term 

C m = e ij e^Mf f w{ a w^ + c.c. (2.6) 



2 Earlier, only the case iV c > 3 when the fermions belong to the essentially complex representation has 
been considered, but the theorem can be easily generalized to the case when the representation is real or 
pseudoreal M). 



(where Ai/p is an arbitrary antisymmetric complex matrix 2Nf x 2Nf) breaks the 
SU(2Nf) symmetry explicitly and gives masses to the Goldstone bosons which are, how- 
ever, small if the matrix elements A4ff> are much smaller than A® CD which we will always 
assume. 

Let us discuss now the symmetry breaking pattern in the case when the fermions 
belong to the adjoint representation of the color group which is real 0, 14, |J. The 
fermion part of the lagrangian^]. is 

£ferm = l $> W, (2.7) 

/=1 



where Nf is now the number of massless Weyl adjoint fermions w^ a The lagrangian 
possesses a U(Nf) symmetry of which the U(l) part is anomalous. In this case the Vafa- 
Witten theorem dictates [13] the condensate to be flavor-symmetric (see also [[F]), 

E 
2 



e^w a Jwf) n = tyf, (2.8) 



(This form is fixed up to an arbitrary unitary transformation from SU(Nf) : 8^ — > 
{UU T yf '.) The formation of the condensate breaks down the symmetry to SO(Nf) 
(Only orthogonal transformations leave the form e a ^w a Jw a ^ invariant). This leads to the 
appearance of 

N 2_ 1 _ NfjNf - 1) _ N f (N f + 1) _ 
f 2 2 

Goldstone bosons. They acquire small masses after switching on the quark mass term 

C« = e^M ff w a Jwf + c.c, (2.9) 

(where A4fp is now a symmetric complex Nf x Nf) matrix. 

Before proceeding further, let us discuss the simplest case Nf = 1 and remind simul- 
taneously some relevant arguments from Q where the reader can find the details. It is 
well known that, in standard QCD with only one light quark, no spontaneous chiral sym- 
metry breaking occurs at all. After explicit breaking of Ua(1) due to the axial anomaly, 



3 The lagrangian (2.7) can also be expressed in terms of the Majorana 4-component fields 



w 
-a 2 w 



(see e.g. |lf 



the symmetry of the quantum lagrangian is just Uy (1) which is not broken spontaneously 
by the formation of a condensate. 

Also in the other theories discussed above there is no spontaneous symmetry breaking 
for one flavor. For the adjoint Nf = 1 case (for zero fermion mass , this theory is just 
N = 1 supersymmetric Yang-Mills), no symmetry is left right from the beginning. For 
N c = 2, Nf = 1, the chiral symmetry is SU(2) but it is not broken by the formation of 
the condensate < e^e af3 Wi a Wjf3 > which is the group invariant. 

That means that the spectrum of all Nf = 1 theories involves a gap and, when the size 
of the system L is much larger than the characteristic scale Aq^ d , the extensive property 
for the partition function holds: 

Z(m, 9) = exp(— e mc (m, 9)V) ~ exp(Smcos^y), (2-10) 

where m <C ^qcd is the quark mass and 9 is the vacuum angle. The normalization 
e TOC (0, 9) = is chosen. The particular combination m cos 9 appears in the first term of the 
Taylor expansion of e vac in m due to Ward identities which enforce all physical quantities to 
depend on the parameters m and 9 in the combination me %e and the requirement of reality 
of e vac . The parameter £ has the meaning of the quark condensate. For defmiteness, we 
have written ( |2.10| ) for theories with fundamental fermions. For adjoint fermions, the 
combination me ie l Nc enters and one should substitute 9 — > 9/N c in ( |2.1(J| ). 

The partition function in the sector with a given topological charge v is given by the 
Fourier integral 

1 r** 
2n Jo 

where I u (x) = I- V (x) is the exponentially rising modified Bessel function. In the adjoint 
case, the integral over 9 extends up to 2ttN c and the same result holds but with rescaling 
v — > v = vN c . The admissible topological numbers v are then integer multiples of 1/N C 
T7|j . The sum rules ( |1 . 1| ) , (|1.4|), and ( |1.6| ) are obtained from the expansion of Z u (m) 



ZJm) = — e- iu8 d9Z(m,9) = UmEV), (2.11) 
2tt Jo 



in m and comparing it with the expansion of (the Euclidean gamma matrices are anti- 
hermitean) 

" (n K)> ■ 

They all coincide for Nf = 1. 
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3 N c = N f = 2. 

When Nf > 2, spontaneous chiral symmetry breaking occurs and Goldstone modes ap- 
pear. The presence of quasi-massless particles in the spectrum brings about long-distance 
correlations which do not allow one to write down the finite volume partition function in 
the extensive form ( |2.10| ). However, the properties of Goldstone bosons are known, and 



one can take them into account explicitly in the path integral for the partition function. 
When L is much less than the inverse Goldstone mass ~ (itiAqcd)" 1 ^ 2 , it suffices to 
take into account only the zero spatial Fourier harmonics of the Goldstone fields and the 
path integral is reduced to a finite dimensional integral over the coset space where the 
Goldstone modes live. 

The case Nf = 2 is particularly simple because of the simple structure of these cosets. 
The integrals can be done explicitly. For N c = 3 with fundamental fermions, the Goldstone 
modes live on SU{2) = S 3 . In the adjoint case, the coset is SU{2) / SO{2) = S 2 . These 
two cases have been treated in ||]. The analysis of the N c = Nf = 2 theory is performed 
along the same lines. 

Again, the coset is simple: 

K = SU{4)/Sp{4) = SO{6)/SO{5) = S 5 . (3.1) 

The corresponding low-energy chiral effective lagrangian belongs to the class considered 
in |18| . It can be formulated in terms of the 6-dimensional vector n with the constraint 



n 2 = 1. When L <C (mAgc'D) -1 ^ 2 , the kinetic term ~ (<9 M n) 2 in the effective lagrangian 
can be neglected, and only the mass term is relevant. If all quark masses are equal and 
real, the mass term is just 2mcos|Sn 6 . As earlier, the particular dependence on the 
vacuum angle 9 follows from the Ward identities. The partition function can be written 

as 

Z(m,6) oc J dn b e 2mcos ^ Vcos ^, (3.2) 

where $ is the polar angle on S 5 (so that dfl^, oc sin $). The integral is elementary and 
we get 



2/ 2 (2mc OB f£VQ 
Z( ' M) = (mcosfSIT ' (X3) 



(we choose the normalization Z(0,9) = 1). The partition function in a given topological 
sector v is the coefficient of e lv0 in the Fourier expansion of Z(m, 9): 



/X\- U ' l + " ] 



[2(|i/ 


+ k))\ 


k\{\u 


+ k)\{\u 


+ k + 2)\(2\v 


+ k)\ 



(3.4) 



with x = mYM . Comparing the coefficients of and x'"' +2 , we arrive at the sum rule 



Q1J) with JVy = 2. 

A more explicit expression for the partition function Z v is obtained is we integrate 
over 9 first. As a result we find 



8 r 

Z v = —\ d$sin 4 $/ 2i ,(2a;cos$). 

37T JO 



(3.5) 



By writing sin $ = |(cos4$ — 4 cos 2$ + 3), this integral can be expressed into known 



integrals f24" |. The result is 

1 



3I u (x) - AI v _ x {x)I v+ i{x) + I v - 2 (x)I v+2 (x) 



(3.6) 



In section 5 we will show that this formula can be generalized to arbitrary Nf in the form 
of the Pfaffian of an anti-symmetric matrix to be defined later. 



4 Sum rules for any Nf. 

Let us consider first the case of adjoint fermions with arbitrary number of flavors iVy The 
coset SU(Nf) / SO(Nf) is spanned by the symmetric unimodular matrices S which can 
be presented in the form S = UU T with U G SU(Nf) . The partition function Z(A4,9) 
represents the invariant integral over the coset which, for convenience, will be traded off 
for the integral over the full group SU(Nf): 

Z(M,9) = J su ^^dUexp{VZRe[Tr{Me i9/NcN fUU T )}}. (4.1) 

Here, dll is the Haar measure on the group, and due to Ward identities, the mass matrix 
M. and the vacuum angle 9 enter in a particular combination M. exp(i9 / 'N C N f) only 
|4|]. The integral ( |4.1|) involves also dummy variables spanning the stability subgroup of 
orthogonal SO(Nf) matrices on which the integrand does not depend. 
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The partition function Z U (A4) in the sector with a given topological charge v is given 
by the Fourier integral 

Z V {M) = — - / e- iv6 d9Z(M,9). (4.2) 
Ztt1\ c Jo 

It can be written as the integral over the full U(Nf) group as [y = vN c ) 

ZJX) = f dU (det U)~ 2D exp (Re(TiXUU T )) , (4.3) 

JU(Nf) v ' 

where X = VEA4. The mass matrix J\A is symmetric and complex. 
The partition function ( |4.3|) is an invariant function of X: 

Z P (V T XV) = (det Vf~ u Z D {X), (4.4) 

where V is a unitary matrix. This enables us to expand the partition function for v > 

as 

Z P (X) = N p [det(X)f(l + aptrX+X + 0(X 4 )). (4.5) 

The normalization factor is denoted by N„. (For v < the first factor in the expansion 
should be replaced by [det(X"'')]~ 1 '.) Such integrals have been considered before in [[| for 
the SUL{Nf)®SUn(Nf) — > SUy(Nf) chiral symmetry breaking scheme. A similar analysis 
will be applied to the present partition function. The matrix X can be decomposed as 

M s 

X = Y,x a t a , (4.6) 

a=l 

where t a is an orthonormal set of symmetric real matrices which are the generators for 
the symmetric unitary matrices. They satisfy the identities 

Tr{t a t b } = ^5 ab , a = l,---,M s , 

M s 

y t a t a = — -, (4.7) 

The total number of generators M s = Nf(Nf + l)/2. The coefficients x a in ( |4.6| ) are 
complex numbers. For any symmetric matrices A and B, we have the identity 

Tr(t a A)Tr(t a B) = -Tr(AB). (4.8) 

0=1 ^ 
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This enables us to show that the partition function (|4.3| ) satisfies the differential equation 

Y,d Xa d Sa Z p = -fz p . (4.9) 



Substituting here the expansion ( |4.5|) , the coefficients a v can be easily found: 

1 



(4.10) 



" 4(|P| + (A> + l)/2)- 

Comparing ( (4.5|) with the mass expansion of the averaged Dirac determinant J3J, one 
immediately gets to the spectral sum rule ( |1.6| ) which agrees with the result obtained by 
using random matrix theory |J. 

The calculation of the sum rule for fermions in the fundamental SU (2)— color represen- 
tation is analogous. In this case the partition function corresponds the chiral symmetry 
breaking scheme SU(2Nf) — > Sp(2Nf). The coset is spanned by the antisymmetric uni- 
modular matrices A which can be presented in the form A = UIU T where U G SU(Nf) 
and / is the symplectic matrix defined in (|2.5| ). In this case the mass enters only in the 
combination M. exp(i6/Nf). The partition function in the sector with a given topological 
charge v is given by the integral over the coset and vacuum angle 9 which can be expressed 
in terms of the integral over the full U(2Nf) group : 

ZJX) = I dU (det U)~ u exp ( -Re(TrXUIU T )) . (4.11) 

Ju(2N f ) \2 / 

The matrix X = VY>Ai is now an anti-symmetric complex matrix. This partition function 
also satisfies an invariance relation 

Z U (V T XV) = (det V) V Z V {X), (4.12) 

(cf. (|4.4|) ), and therefore can be expanded in group invariants. The expansion is a little 
bit modified as compared to fl4.5|) and for v > it has the form 



Z U (X) = N u [Pi(X)] u (l + ^a.tiX^X + 0{X i )). (4.13) 



where Pf(X) = ydet(X) is the Pfaffian of the antisymmetric matrix X. (For v < the 
first factor in this expansion has to be replaced by [Pf (X^)]~" .) 

The matrices UIU T are anti-symmetric complex unitary matrices. Their generators 
are anti-symmetric real 2Nf x 2Nf matrices which we use as a basis for the expansion of 
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X: 

M aa 

X=Y,Xat a . (4.14) 

a=l 

In this case the total number of generators is M as = Nf(2Nf — 1). Instead of ( |4.7| ) we 
have the relations 

Tr{t a t b } = h ab , a = l,---,M as , 

EM as 
t a t a = —, (4.15) 

and ( [4.8| ) is also valid for anti-symmetric matrices A and B. In this case the sum over 
a runs over the anti-symmetric generators. The partition function ( [4. 1 1| ) satisfies the 
differential equation 

Nf 

Y^d Xa d Sa Z y = -lz v . (4.16) 
(cf. (|4.9| )). Substitution of the expansion of Z u in powers of X yields 

a » = 4(1,1+2^-1) - (417) 

The corresponding sum-rule has the form ( |1.4| ) which agrees with the result obtained from 
random matrix theory. 

The sum rules ( |1 . 1| ) , ( |1.4| ), ( |1.6| ) can be written universally as 

\&oW v 4[|H + (dim(coset) + l)/iV / ]' l ' ] 

with the rescaling v — > v and counting in the sum only one eigenvalue of each degenerate 
pair in the adjoint case. 



5 Partition function for equal masses 

In this section, we evaluate the partition function for the special case when all masses are 
equal to m. We use the notation x = mVE. First, we consider the case of adjoint fermions. 
The matrix UU T in (4.1) is a symmetric unitary matrix which can be diagonalized by 
an orthogonal matrix. To evaluate the integral we use the eigenvalues exp(z#fc) and the 
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eigenvectors of UU T as new integration variables. The Jacobian of this transformation 
has been used extensively in random matrix theory. It has the form 



J s (e 1 ,---,9 Nf )oc 1 [[\e^-e iei \. (5.1) 



k<i 

The partition function can be rewritten as 

N 



Zv{x) - ^ 



r Nf/2 f7T f7T m / \ 



(5.2) 



The normalization is such that Z(0) = Zq(0) = 1. Integrals of these type have been 
studied extensively by Mehta pOj, |2l| . The result can be expressed as a Pfaffian (see 
appendix A) 

ZJx) = f —— Pi(A). (5.3) 

uK ' 2 N fT(l + N f /2) K ' K ' 

For even Nf, the matrix elements of the anti-symmetric matrix A are given by Q 

— / — e(6 — <p) exp(i(p<j) + q6)) exp(x cos0 + x cos 9 + iv((j) + 9)). (5.4) 
-7t 2n J -it 2ix 

The indices p and q run between — + | and ^ — \. 

For v = and x — 0, the Pfaffian in (5.3) can be calculated readily. The result is 

2 N f /2 

Pi ( A )l=0,x=0 = n N f /2( Nf _ -QH' ( 5 - 5 ) 

which gives the correct normalization for the partition function. 

The simplest non-trivial case is Nf = 2 which describes two Major ana adjoint flavors 
(or one Dirac flavor). In this case, one can derive 

7T r 2 l p l / r 2 \ 

ZJx) = -A 11= — 1 + —f + 0(x A ) , (5.6) 

v ' 2 ~ 2 (2 i7 + 1 ! ^ 2z? +3 V 7 ' K ' 



4 



Expanding e(0 — <fr) in the Fourier series on the interval — 2ir < 9 — <fr < 2tt, one can present the 



integral (5^4) in the form 

1 °° 1 

= ~ X! lTT 1 I v+p+k+±{ X ) I v+q-k-} 1 { X )- 



71 u k 1 2 
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which agrees with results obtained in ||. For odd values of Nf the partition function is 
equal to the Pfaffian of the matrix A supplemented by an additional row and column (see 
appendix A) 

7T N f /2 N f \ (A c\ 

Zp = 2 N fT(l + N f /2) Pi \-c T ) ' (5 ' 7) 
where c is a vector of length Nf defined by 

— exp(ip8) exp(x cos 6 + iv6). (5.8) 

Again, the Pfaffian can be calculated easily for v — and x = 0: 

2 (7V/-l)/2 

Fi ( A )\»=o,x=o = n (N f -i)/2( Nf _ iyv ^ 5 " 9 ) 

which results in a partition function that is normalized to one. 

The small x expansion of the partition functions ( |5.3| ), ( p.7| ) starts from the term 
oc x Nf ^ which corresponds to the presence of Nf\i>\ pairs of fermion zero modes in the 
gauge field background with topological charge v. For x ^> 1, 

Z p oc exp(N f x) (5.10) 

for all V which means that, in the infinite volume limit, the contributions of all topological 
charges are of the same order |3J]. 

For N c = 2 with fundamental fermions, we proceed in exactly the same way. In this 
case the anti-symmetric unitary matrix UIU T can also be brought into the canonical form 
by an orthogonal transformation. The result is an anti-symmetric tri-diagonal matrix with 
off-diagonal matrix elements ±exp(z#fc). Again we use the angles 9k as new integration 
variables. The Jacobian of this transformation is 

^(^•••.^ocnie^-e^l 4 , (5.11) 

k<i 

leading to the partition function 

z - - m -w L ■ ■ ■ L S * n - ^ i 4 ^ (» ? - «o + * e *>) . 



fe=i fc</ \ & 



(5.12) 
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The normalization factors have been chosen such that Zq(0) = 1. Also this partition 
function can be rewritten as a Pfaffian (see appendix B): 

1 



" (2JV>-1)!! 
where the matrix elements of A are given by 



Pf(A), 



(5.13) 



/7T dt> 
— exp(i(q + p)9) exp(x cos(6 l ) + iu6 
-■K 2ir 



(q - p)Ip +q+u , 



(5.14) 



where the indices p and q run between —(Nf — \ ) and Nf — |. Note that A is a 2Nf x 2iV/ 
dimensional matrix. 

For x < 1, Z v (x) oc a?^'"! which agrees with the zero modes counting. The large x 
asymptotics of Z u (x) is given by ( |5.10| ) and is universal for all v. 

As an example, let us consider Nf = 2. In this case the matrix A is given by 



A 



I I v _ 2 (aO 2J I/ _i(x) 3I u {x) \ 

-h- 2 (x) 4(x) 2I v+1 (x) 

-2I v ^{x) -I v (x) J, +2 (x) 

V -3I v (x) -2I u+1 (x) -I u+2 (x) 



(5.15) 



The Pfaffian of this matrix is given by 



Pf(A) = 3I 2 u (x) - U v _ x {x)I v+1 {x) + I u . 2 (x)I u+2 (x), (5.16) 

and Z v {x) coincides with our previous result ( |3.6j ) derived by performing the integral over 
the angles directly in ( |3.5| ). Expansion of Z v {x) in powers of x reproduces the sum-rule 
derived in previous sections. 



6 Discussion 



Depending on whether the representation of the color group is real (SU(N C ), N c > 2 
with adjoint fermions) , complex (SU(N C ), N c > 3) or pseudoreal (SU(2) with fun- 
damental fermions), chiral symmetry breaking is realized in a different way leading to 
different Goldstone sectors. They are parameterized by SU(Nf)/SO(2Nf), SUn(Nf) x 
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SUL{Nf) I SU(Nf) and SU(2Nf)/Sp(2Nf), respectively. For each of the three cases we 
have determined explicitly the dependence of the partition function in finite space-time 
volume on the fermion mass m. As was shown in Q , the comparison of the chiral ex- 
pansion of these partition functions with the chiral expansion of the Dirac determinant 
in the original path integral leads to sum rules for the inverse powers of the eigenvalues 
of the massless Dirac operator in QCD. In other words, we have very specific correlations 
between the eigenvalues of the Dirac operator. 

The same triality is known from random matrix theory. The matrix elements of the 
three classical random matrix ensembles are real, complex or quaternion real. These 
ensembles can be generalized to include the chiral structure of QCD, and it is possible to 
derive the spectral density and its correlation functions. Previously it had been verified 
that for complex representations of the color group, this leads to exactly the same sum 
rules as obtained from the finite volume partition function [2ql- In the framework of 



random matrix theory it is straightforward to derive sum rules for the two other ensembles 
as well ||. (The theories with real, complex , and pseudoreal fermions correspond to 
setting (3 = 4, (3 = 2, and (3 — 1 in the universal formula (14) of Ref. 0). 

This leads to the question whether the sum rules for the two remaining cases are 
indeed the ones that follow from the effective partition function that is based on general 
assumptions only. This question has been answered affirmatively in this paper. Since both 
approaches are based on symmetry this result did not come as a surprise. Apart from 
this, several other new results were derived. In particular, we found explicit expressions 
for the partition function for equal quark masses which generalize results already known 
for the complex case. This does not yet answer whether all spectral correlations obtained 
from random matrix theory are also properties of the QCD partition function. However, 
from the study of spectra of classically chaotic systems, we know that correlations on the 



order of one or several level spacings are universal fll9| . It are precisely such correlations 
that are responsible for the spectral sum rules. This leads to the conjecture that the 
microscopic spectral density is universal. 

In this paper, we always assumed that the parameters of the theory were in the range 

Aqcd « L « m- 1 ~ {mK QCD y 1 ' 2 . (6.1) 
17 



The upper bound for the length of the box L was essential because, if L would be larger 
than the Compton wavelength of the Goldstone modes, ~ (mAg^) -1 / 2 , nonstatic modes 
of the Goldstone fields would be excited and the path integral could not be reduced to a 
finite dimensional group integral. The calculation is also possible for larger L by exploiting 
the fact that the characteristic fluctuations of Goldstone fields are small in that case |26[ 



But, as we are mainly interested in the derivation of sum rules ( |1.4j) and ( |1.6| ) which 
are formulated for theories with massless fermions, results in the range ( |6.1| ) are quite 
sufficient for our purposes. 

From the practical point of view, the results for the theory with N c = Nf = 2 analyzed 
in section 3 are perhaps the most important. This theory bears all features of standard 
QCD but involves fewer degrees of freedom so that it will be easier to confront the sum 
rules with lattice simulations. 

In this paper, the effective partition functions were postulated on the basis of chiral 
symmetry and other general arguments. Our results suggest that the static limit of the 
finite volume partition function can also be derived directly from random matrix theory. 
Work in this direction is under way. 
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Appendix A 

In this appendix we calculate integrals of the type 

r w r-K N Af), N 

P (u)= ... n«(^ni^-^i. (A.i) 

J - w J - W k=l An k<l 

where u(8) is an arbitrary function of 9. We only give the main steps and refer to Dyson's 
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original paper [22 or the book of Mehta pOf for further details. Because the integrand is 
a symmetric function of the integration variables, the integration region can be restricted 
to the domain D defined by — 7r <#i<---<#yv<7rat the expense of a factor N\. For 
angles inside this domain we have the identity 



n 

k<l 



JV(JV-l) 



— det{0 p (^)}, 



(A.2) 



where the determinant is over anJVxiV matrix with matrix elements given by 

N-l N-3 N - 1 



2 ' ' 2 



(A.3) 



To proceed we use the method of integration over alternate variables originally due to 
|25|j . The result for even N is 

S i-iv (P 1 ?.) 4> i-n (P 1 ?.) •■• S i-N (9iv) 4> i-n (Oiv) 

p(u) = i —N\ / J] ——u(9k) 
Jd, 2tt 

k even 



S N-l (9<).) (f) N-l {Ol) •"■ S N-l {On) 4> N-1 (Ojv) 



(A.4) 



The quantities S p (9k) are defined as 



S, 



p\Vk> 



u((f))(j) p ((t)) 



2vr' 



(A.5) 



For odd N, the last column is Si-n (it), ■ ■ ■ , Sn-i (it). The integrand in (A.4) is a sym- 
metric function of the integration variables, which permits us to extend the integration 
range of each variable to [— n, tt] and use the integration theorem of appendix (A. 7) of ref. 



20 . As a result we obtain 



p{u) = N\ Pi{A), 



(A.6) 



where the matrix elements of the anti-symmetric matrix A are given by 

<•* de, 



/ n (16 
-[S p (6)<i> q (6) - S q 
-7r 2n 



2 L L ~ ^X^M^M - M<P)Md)l (A.7) 
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(the identity i ( 2 ' = (— 1}1 ( even AT, has been used). The Pfaffian of an N x N (N 
even) anti-symmetric matrix is defined as 

Pf ^ = (N/2)\ ^ si g n ( (T ) A -(i)-(2) ■ ■ ■ Mn-iMN) , (A.8) 
where the sum is over all permutations a of {1,2, ■ ■ -N}. For Pfaffians we also have a 



Laplace expansion (see p. 393 of p7| ) 



Pf(A) = (-l) r+1 [A rl Pf(A lr , lr ) - A r2 Pf(A 2r , 2r ) - • • ■ + (-l) N+1 A rN Pf(A Nr , Nr )], (A.9) 

where the matrix A^kr is obtained from matrix A by deleting the both the fc'th and the 
r'th rows and columns. 

For odd N, we expand the determinant with respect to the last column. Then the 
theorem from appendix (A. 7) of [^(J can be applied to each of the minors after extending 



the integration domain of each of the integration variables to [— 71", tt]. The extra factor 
1/[(N — l)/2]! cancels against the same factor from the integration theorem in ||20|| . Using 



the inverse of the Laplace expansion (|A.9|) we find (with the use of the identity i ( 2 ' = 
(i)V, odd N) 

M = N\Pi(_ A qT (A. 10) 
where the matrix elements of A defined above and 

N—l 

c p = (-1) — S p (n). (A.ll) 



Appendix B 

In this appendix we calculate the integral 

P (u)= ... n^)|rni^-e^| 4 . (B.l) 

J - w J - n k Z7r k<i 

According to Mehta, the Jacobian can be rewritten as a determinant of a 2N x 2A^ matrix: 

N 

II V ldk - = detixpiejlpXpiOj)}, (B.2) 

k<l 
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where 

toft) = p = I. (B.3) 

In (B.2), the rows of the determinant are indexed by p, whereas the columns are indexed 
by j. According to the theorem in appendix (A. 7) of [2(| the integral (B.l) can now be 
expressed Pfaffian 

p(u) = iV!Pf(A), (B.4) 
where the matrix elements of A are given by 

— (x P (0fc ? (0)] - x g (0)bx#)D • (B.5) 

-7T 27T 



21 



References 

[1] Reprint volume 'Vacuum Structure and QCD Sum Rules' (M. A. Shifman, ed.), North- 
Holland, 1992. 

[2] Lattice 1992: Proceedings of the international symposion on lattice field theory, (J. 
Smit and P. van Baal eds.), North-Holland, Amsterdam 1993. 

[3] D.I. Diakonov and V.Yu. Petrov, Nucl. Phys. B245 (1984) 259; E. Shuryak and J. 
Verbaarschot, Nucl. Phys. B341 (1990) 1; Nucl. Phys. B 410 (1993) 55. 

[4] H. Leutwyler and A. Smilga, Phys. Rev. D46 (1992) 5607. 

[5] E. Shuryak and J. Verbaarschot, Nucl. Phys. A560 (1993) (306). 

[6] J. Verbaarschot, The spectrum of the Dirac operator and chiral random matrix theory: 
the threefold way, Phys. Rev. Lett. (1994) (in press). 

[7] M. Peskin, Nucl. Phys. B175 (1980) 197. 

[8] C. Vafa and E. Witten, Nucl. Phys. B 234 (1984) 173. 

[9] D. Weingarten, Phys. Rev. Lett. 51 (1983) 1830. 
[10] E. Witten, Phys. Rev. Lett. 51 (1983) 2351. 
[11] S. Nussinov and M. Spiegelglas, Phys. Lett. B181 (1986) 134. 
[12] S. Coleman and E. Witten, Phys. Rev. Lett. 45 (1980) 100. 
[13] M. Shifman, Private communication. 

[14] S. Dimopoulos, Nucl. Phys. B168 (1980) 69; M. Vysotskii, I. Kogan and M. Shifman, 
Sov. J. Nucl. Phys. 42 (1985) 318. 

[15] D.I. Diakonov and V.Yu. Petrov, in 'Quark Cluster Dynamics' , Proceeding of the 
99th WE-Heraeus Seminar, Bad Honnef, 1992 (eds. K. Goeke et al.), Springer 1993. 

[16] P. Ramond, Field Theory: A Modern Primer, Benjamin/Cunnings, Reading, MA, 
1981. 

[17] E. Cohen and C. Gomez, Phys. Rev. Lett. 52 (1984) 237; A. Zhitnitsky, Nucl. Phys. 
B340 (1990) 56. 

[18] P. Hasenfratz and H. Leutwyler, Nucl. Phys. B343 (1990) 241. 

[19] O. Bohigas, M. Giannoni, in 'Mathematical and computational methods in nuclear 
physics', J.S. Dehesa et al. (eds.), Lecture notes in Physics 209, Springer Verlag 
1984, p. 1; O. Bohigas, M. Giannoni, and C. Schmit, Phys.Rev.Lett. 52, 1 (1984); 
T. Seligman, J. Verbaarschot, and M. Zirnbauer, Phys. Rev. Lett. 53, 215 (1984); 



22 



T. Seligman and J. Verbaarschot, Phys. Lett. 108A (1985) 183; B.Simons, A. Szafer 
and B. Altschuler, Universality in quantum chaotic spectra, MIT preprint (1993); E. 
Brezin and A. Zee, Nucl. Phys. B402 (1993) 613; H. Weidenmiiller, in Proceedings 
of T. Ericson's 60th birthday. 

[20] M. Mehta, Random Matrices, Academic Press, New York, 1967. 

[21] M. Mehta, Random Matrices, Academic Press, San Diego, 1991. 

[22] F.J. Dyson, J. Math. Phys. 3 (1962) 140. 

[23] J. Verbaarschot and I. Zahed, Phys. Rev. Lett. 70 (1993) 3852. 

[24] I. Gradshteyn and I. Ryzhik, Tables of Series, Integrals and Products, Corrected and 
enlarged edition (Academic, San Diego, 1980). 

[25] M. Mehta, Nucl. Phys. 18 (1960) 395. 

[26] J.Gasser and H.Leutwyler, Phys. Lett. B 184 (1987) 83; 188 (1987) 477. 

[27] T. Muir, A treatise on the theory of determinants, 1933 (reprinted by Dover, New 
York, 1960). 



23 



